Multicritical points for the spin glass models on hierarchical lattices 
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The locations of multicritical points on many hierarchical lattices are numerically investigated by 
the renormalization group analysis. The results are compared with an analytical conjecture derived 
by using the duality, the gauge symmetry and the replica method. We find that the conjecture does 
not give the exact answer but leads to locations slightly away from the numerically reliable data. We 
propose an improved conjecture to give more precise predictions of the multicritical points than the 
conventional one. This improvement is inspired by a new point of view coming from renormalization 
group and succeeds in deriving very consistent answers with many numerical data. 
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I. INTRODUCTION 

The problem of spin glass, which is one of the most 
challenging subjects in statistical physics, has been an- 
alyzed extensively by the mean field theory but has 
not been sufficiently understood for finite dimensional 
systems. P, [2, S| Most approaches to the difficult prob- 
lem of finite dimensional spin glasses rely on approximate 
techniques, numerical simulations and phenomenological 
theories. Q 

A lot of important facts have been established by such 
approaches to finite dimensional spin glasses. Neverthe- 
less, it is very important to derive exact or rigorous re- 
sults to check validity of approximate approaches. A use- 
ful method along this line is the gauge theory, which en- 
ables us to find a special subspace of a phase diagram 
for spin glass models known as the Nishimori line. [J, Q 
In this subspace, we can calculate the exact value of the 
internal energy and evaluate the upper bound of the spe- 
cific heat. It is shown rigorously that the Nishimori line 
runs through the ferromagnetic phase and the param- 
agnetic phase. Moreover it is expected that the multi- 
critical point, where the phase boundaries between the 
spin glass phase, paramagnetic phase, and ferromagnetic 
phase merge, is located on the Nishimori line.[6| 

One of the recent developments using the gauge theory 
is a conjecture of the exact location of the multicritical 
point for spin glasses, especially in two dimensions. 0, Q 
The prediction is very close to numerical results, @ and is 
considered to be very useful in the analysis of numerical 
data for critical exponents. 

The essential part to derive the conjecture consists of 
duality and the replica method as explained below. Du- 
ality is a useful tool to obtain the exact location of the 
transition point for spin systems without disorder. Be- 
cause the spin glass models have quenched disorder, we 
cannot directly apply the duality to spin glass models. 
Nevertheless, by the replica method, the problem reduces 
to non-random systems to which we can apply duality. 
We then assume that a single relation gives the multi- 
critical point similarly to the case of the duality relation 



on the transition point for the pure Ising model in two 
dimensions. This is one of the hypotheses on the con- 
jecture. In addition we expect that the above-mentioned 
relation for the multicritical point is satisfied even when 
the replica number goes to zero to analyze systems with 
quenched randomness. 

The validity of assumptions to derive the conjectures 
as mentioned above has not been established rigorously. 
Nevertheless the conjecture has given predictions very 
close to many independent numerical results, ffi, M. H, flO. 
1 1 II ] However, Hinczewski and Berker found several exam- 
ples by the exact renormalization analysis in which, es- 
pecially for the hierarchical lattices, the conjecture did 
not give good predictions. It is expected generally 
that the renormalization group analysis on hierarchical 
lattices gives exact results. Therefore such discrepancies 
found by the renormalization group analysis on hierar- 
chical lattice should be taken seriously for the conjecture 
even though the amount of discrepancies is small. If these 
discrepancies are genuine, we have to consider the reason 
why there are cases for which the conjecture does not 
work well. Conversely, we would like to know why the 
conjecture always gives accurate results if not exact. It is 
also desirable to improve the conjecture to predict more 
precise location of the multicritical point. 

These observations have given us motivations to in- 
vestigate many more hierarchical lattices to check if the 
locations of the multicritical points are away from pre- 
dictions by the conjecture. We employ the technique by 
Nobre[l3| to examine the phase transitions and to esti- 
mate the locations of the multicritical points on several 
hierarchical lattices. Then, we discuss reasons why there 
are such discrepancies for the cases of hierarchical lat- 
tices and propose a method to improve the conjecture to 
derive more precise locations of multicritical points. 

The presented paper is organized as follows. In Sec. 

II, we review previous results for the location of the mul- 
ticritical point, and show some examples of the discrep- 
ancies between the conjecture and the numerical results. 
We explain the properties of hierarchical lattices in Sec. 

III. In Sec. IV, we explain Nobre's method to examine 
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the phase transition and then estimate the values of the 
locations of the multicritical points for several cases. We 
find some discrepancies between the conjecture and the 
numerically estimated results here. Therefore we have 
to consider the reason why there are some cases that 
the conjecture does not work well. The conjecture relies 
on the replica method, and we assume that the analyt- 
ical continuation of the replica number to zero does not 
cause troubles as in most cases studied so far. Therefore 
we investigate the replicated systems for the ±J Ising 
model in Sec. V. Then we consider improvement of the 
conjecture in Sec. VI and show successful results to pre- 
dict more precise locations of the multicritical points in 
Sec. VII. In Sec. VIII, discussions and future outlook 
are given. 



Type Conjecture Numerical result 



SQ ±J Pc = 


0.889972[7, 8] 


0.8900(5)[14] 






0.8894(9) [15] 






0.8907(2) [16] 






0.8906(2) [17] 






0.8905(5) [18] 






0.8907(4) [19] 


SQ Gaussian Jo/ J 2 


= 1.021770[7, 8] 


1.02098(4) [19] 


TR ±J Pc = 


= 0.835806[11] 


0.8355(5) [14] 


HEX ±J Pc = 


= 0.932704[11] 


0.9325(5) [14] 



TABLE I: Comparisons between the conjectured values and 
the numerical results. SQ denotes the square lattice, TR 
means the triangular lattice, and HEX expresses the hexago- 
nal lattice. 



II. MODEL AND CONJECTURE 

We study the random-bond Ising model, defined by the 
Hamiltonian. 



H 



(ij) 



(1) 



where <7j is the Ising spin taking values ±1, and Jy de- 
notes the quenched random coupling. In this paper, we 
consider two types of distribution functions for Jy, the 
± J model and the Gaussian model. 

For this random-bond Ising model on two-dimensional 
lattices, a method to predict the precise location of mul- 
ticritical point has been proposed. 0, H, H, OH [H| This 
method relies on the duality and the replica method ap- 
plied to spin glass models with gauge symmetry. It has 
been considered to be a conjecture for the exact location 
of the multicritical point for systems satisfying certain 
conditions like self duality. According to this conjecture, 
the exact location of the multicritical point for the ± J 
Ising model on the square lattice is determined by a single 
equation as follows, 

-plog 2 p- (1 -p)log 2 (l -p) = ~, (2) 

where p is the probability of Jy = J > for the ± J Ising 
model. The left-hand side of this equation is the binary 
entropy and will be written as H(p). We obtain the value 
p c = 0.889972 (ss 0.8900), solving this equation. This 
result is in reasonable agreement with existing numerical 
results as shown in Table fl] It is also possible to obtain 
the location of the multicritical point for the Gaussian 
model with the average Jo and the variance J 2 from the 
following equation, 0,1a] 



dJijPiJij) log 2 {1 + exp(-2/3J y )} = i (3) 



where P(Jy) is the Gaussian distribution function, and 
(3 satisfies the condition of the Nishimori line (3 = Jo/ J 2 . 



We will write the left-hand side of this equation as 
Hg(Jq/J 2 )- The solution of Eq. J3]) for the Gaussian 
model is Jo/ J 2 — 1.021770, which is also compared with 
the existing numerical result in Table [IJ It is not easy to 
determine from these data whether the conjecture actu- 
ally gives the exact result. 

Equation ([2]) applies also to models defined on other 
self-dual lattices. The phase diagram of a self-dual hi- 
erarchical lattice has been numerically investigated by 
Nobre.[l3| According to his result, the multicritical point 
is located near the conjectured value, p c — 0.8902(4). 

In addition, the conjecture also works on mutually dual 
pairs of lattices. [13] In this case, we obtain the relation- 
ship between the locations p\ and P2 of the multicritical 
points for the mutually dual pairs as follows 



Hipx) + H{p 2 ) = 1. 



(4) 



This relation is supported by a consistent result within 
its numerical error bar for the ±J Ising model on the 
triangular and hexagonal lattices as H(pi) + H(p2) — 
1.002 (3), [3] where p\ and pi denote the location of the 
multicritical point on the triangular p\ — 0.835806 and 
the hexagonal P2 = 0.932704 lattices, respectively. [ll| 

However there are cases in which the relation (fj[|) and 
the numerical results for three mutually dual pairs of 
hierarchical lattices show derivations by large amounts 
close to 2%, H( Pl ) + H{p 2 ) = 1.0172, 0.9829, 0.9911. [H 
The technique by Hinczewski and Berker in Ref. [12] is 
based on an exact calculation through the renormaliza- 
tion group analysis on hierarchical lattices. 

These results motivated us to study other hierarchical 
lattices to see if the conjecture gives exact solutions. If 
it does not, the next question is why the prediction of 
the conjecture falls very close to numerical estimates in 
all cases. To verify these points, we evaluated the rela- 
tion ([J| and its Gaussian version for other five mutually 
dual pairs of hierarchical lattices. In addition, we also 
reexamined Eq. ^ for the ± J Ising model and Eq. ([3]) 
for the Gaussian model on several self-dual hierarchical 
lattices including the case investigated by Nobre. 
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III. HIERARCHICAL LATTICE 

In this section, we introduce the hierarchical lattice 
and renormalization group on it. [2(| [U, [22j The renor- 
malization group analysis on the hierarchical lattice is 
an exact technique to obtain the location of the transi- 
tion point, though it is difficult to obtain such an exact 
solution on regular lattices. 

In the present paper, we investigate phase transi- 
tions for three mutually dual pairs of hierarchical lattices 
shown in Fig. Q] studied by Hinczewski and Berker[l2} 
and additional five pairs shown in Fig. O We also ex- 
amine phase transitions in several self-dual hierarchical 
lattices in Fig. [3] The scale factor 6 denotes the length 
of the unit of the hierarchical lattice. We examine phase 
transitions and estimate the location of the multicritical 
point, for 6 = 2,3, 4, 5, and 6 on the self-dual hierarchical 
lattices. 

The renormalization group on these hierarchical lat- 
tices consists of two basic steps, which are known as the 



bond moving and decimation as in Fig. [4j The black bold 
bond denotes the renormalized bond after bond moving, 
and the white bold bond expresses the renormalized bond 
after decimation. Construction of a hierarchical lattice 
starts from a single bond, and we iterate the process to 
substitute the single bond with a unit cell of more com- 
plex structure as in Fig. O 

Because a hierarchical lattice has an iterative structure 
consisting of unit cells as shown in Figs, [lj [2] and [3j we 
again obtain the same structure after we trace out degrees 
of freedom on each unit cell in renormalization group cal- 
culations, which are the inverse processes of the construc- 
tion. Therefore our task is to evaluate recursion relations 
of couplings on bonds, which relate sets of the couplings 
{K^j,} after renormalization with {K^ ^} before renor- 

(r) 

malizations. The superscript of K\-' denotes the step of 
renormalization. For example, the explicit recursion re- 
lations for the 6 = 2 self-dual hierarchical lattice in Fig. 
are 



e 



El zr(r-l) , iv-(r-l) , tHr-l) , tAt-1) . clr-l) 

exp f Kfa 'ai+Kfa V 2 + K{ 2 'o\.<j% + JQ 3 Vi+JQ 3 V 2 

to} 

exp IKfa Vi + K( )2 '<7 2 + K{ 2 '<ti<7 2 - 0i - K 23 ' a 2 



(5) 
(6) 



Here x , which we call the principal Boltzmann factor, 
expresses the local (bond) Boltzmann factor for paral- 
lel spins on the ends of renormalized bonds. The recur- 
sion relations ([5]) and © yield the renormalized principal 
Boltzmann factor and that for antiparallel spins on the 
ends of renormalized bonds, respectively. The summa- 
tion in the exponent is over all interactions in the unit 
cell. The indices of couplings express bonds as labeled in 
Fig. El 

The partition function Z s for a hierarchical lattice after 
s-step construction is generally evaluated as 

= 4 1) " ri, ^- 1 ({^ ) }) 
= 4 2)< \- 2 ({4 2) }) 

= ... = 4 S) B zo({Kl?}), (7) 

where Ng represents the number of bonds at the sth 
step of construction and z s _ r is the partition function 
after r-step renormalization, which is normalized by 

B , namely z s _ r = (^o^) Because of 

this normalization, the value of z s _ r is simply 2 Ns ~ r in 



the high-temperature limit and becomes 2 in the low- 
temperature limit. [9( Here N a — r denotes the number of 
sites after r steps of renormalization for s steps of con- 
struction. In addition, one notices that the number of 
the construction step s decreases effectively at each step 
of renormalization. 

We obtain the free energy per site as, 

-/?/ s tt4 0) }) 

at(0) 1 

= ^iog4 s) ({4 s) » + ^ kg *>({*&>}), 

(8) 

where iV^ = 1. Therefore the free energy per site of a 
model on a hierarchical lattice is generally written as, in 
the thermodynamic limit s — * oo, 

lim -pf s {{K\f}) 

s— >oo J 

= i{^log^({Ai-)}) + ^lo g *({^})}. 

(9) 

The last term in this equation can be calculated for the 
periodic and free boundary conditions by the fact that 



l-(a) l-(b) 2-(a) 2-(b) 3-(a) 3-(b) 



FIG. 1: Three mutually dual pairs of hierarchical lattices investigated in Ref. [12]. The numbers 1, 2, and 3 denote these three 
mutually dual pairs of hierarchical lattices studied in the presented paper. 




o 

4- (a) 4-(b) 5- (a) 5-(b) 




(a) 6-(b) 7- (a) 7-(b) 8- (a) 8-(b) 



FIG. 2: Additional mutually dual pairs of hierarchical lattices studied in the presented paper. The solid lines denote bonds 
replaced by the renormalized bonds at each renormalization. Bonds shown dashed stay unrenormalized. 



6 = 2 



6 = 3 




6—1 bonds 



6—1 bonds 




FIG. 3: Self-dual hierarchical lattices. After renormalization 
of bond moving and decimation, the self-dual hierarchical lat- 
tices become the structure like the Wheatstone bridge. 









o 



o 



6 



8 = 2 



FIG. 5: Construction of one of the self-dual hierarchical lat- 
tices. The number s denotes the construction step. 



FIG. 4: Two renormalization steps. The left-hand side is 
bond moving. The right-hand side is decimation. 



the hierarchical lattice becomes a single bond after suffi- 
cient steps of the renormalization as 



2 (periodic) 
2{l + exp(-2if,J o) ))} (free). 

(10) 



This is negligible due to Noo — > oo for the case of the pe- 
riodic boundary condition. This choice of the boundary 
condition does not affect the results. Therefore only the 
first term log^Q 00 ^ in Eq. §§§ is significant. This quantity 

is a function of of renormalized couplings {.K'jj }, which 
in general obeys a non-trivial distribution in quenched 
random systems. In the next section, we observe the flow 
of the renormalized couplings {ify } using a stochastic 
technique by Nobre[H| to investigate phase transitions 
on the hierarchical lattices, and estimate the location of 
the multicritical point for the ±J Ising model and the 
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Gaussian model. 



IV. QUENCHED SYSTEMS 

In Nobre's implementation of renormalization group 
for disordered systems on hierarchical lattices, [13[ we first 
produce a sample pool of interactions, following the ini- 
tial distribution. For example, our analysis starts from 
preparation of a pool following the ± J or Gaussian dis- 
tribution. Then we randomly choose bonds from this 
sampling pool and form a unit cell of the hierarchical lat- 
tice under consideration. In this unit cell, we carry out 
the renormalization calculation using adequate equations 
such as Eqs. (O and ^ and obtain renormalized inter- 
actions. Iterating this procedure for the other bonds, we 
obtain another pool consisting of the renormalized inter- 
actions, which follows a new type of distribution function 
of the renormalized interactions. Using this renormalized 
distribution, we reproduce a pool of the renormalized in- 
teractions, and iterate the above procedures while ob- 
serving the moments of interactions (Kij) and (Kfj) at 
every step where (• • • ) means the configurational aver- 
age over the renormalized distribution. If (Kij) goes to 
infinity, it is considered that the renormalization flow of 
the sampling pool is attracted toward the ferromagnetic 
fixed point in the interaction space. On the other hand, 
when (Kij) goes to zero, two possible scenarios are con- 
sidered. To distinguish two different scenarios, we have 
to observe (Kfj). If this moment goes to infinity, it is a 
signal that the sampling pool is attracted toward the spin 
glass fixed point. Otherwise, we find a signal that (Kfj) 
falls zero, then the sampling pool goes to the paramag- 
netic fixed point. An additional scenario is seen in the 
present study, which has not been investigated by Nobre's 
method yet. As shown in Fig. [2l the hierarchical lattices 
of type 5, 6 and 8 include a part of interactions following 
the initial distribution function in the unit cell (see the 
dashed lines). These interactions induce the possibility 
of a fixed line like the Kosterlitz-Thouless (KT) phase 
that the sampli ng p ool goes toward neither (Kij) — > oo 
nor (K^) — > 0.[23| This fixed line can be detected by 
\K^-K%-V\^0. 

We introduced 2,000,000 bonds as the set of a sam- 
pling pool, and prepared 1,000 sampling pools in the 
present study, except for the lattice number 5 (a) and 
(b) in Fig. [2l whose sampling pool has 1, 800, 000 bonds. 
We observed the resulting phases after 30 steps of renor- 
malization iterations. For the hierarchical lattices with 
the possibility of a KT transition, we carry out the renor- 
malization of 50 more steps than the other hierarchical 
lattices. Because the investigations are carried out for 
the hierarchical lattices of finite size and with finite num- 
ber of bonds in the sampling pool, we cannot find clear 
boundary expressing the phase transition. In fact, for a 
given lattice, some sampling pools go to the ferromag- 
netic fixed point (or KT phase) and others are attracted 
toward the paramagnetic fixed point (or KT phase). We 
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FIG. 6: Results for the ±J (upper panel) and Gaussian 
(lower panel) Ising models along the Nishimori line on the 
self-dual hierarchical lattice with 6 = 3, see Fig. [3] The 
conjecture states that the multicritical points are located at 
p c = 0.889972 and J /J 2 = 1.02127, respectively. The white 
square denotes the probability of the paramagnetic phase and 
the black one represents that of the ferromagnetic phase. The 
error bars represent l/vT000 reflecting the number of the 
sampling pools. 



obtained the probabilities of appearance of each phase as 
a result by this method. For example, the results for the 
± J Ising model and Gaussian model on Nobre's self-dual 
hierarchical lattice depicted in Fig. [3]are shown in Fig. [6] 
For the other hierarchical lattices, we obtain similar re- 
sults to these plots. We explain the obtained plot in Fig. 
[7] below, which concerns the error bars for these investi- 
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Lattice 



2H(p c 



0.8896 0.8898 0.8900 0.8902 0.8904 0.8906 
P 



FIG. 7: Size effect for the ±J Ising model on the self-dual 
hierarchical lattice with 6 = 3. The black and white marks 
denote the probabilities of the paramagnetic and ferromag- 
netic phases, respectively. The symbols □ and ■ are for 10 3 
sampling pools with 10 6 bonds, A and A are for those with 
5 x 10 6 bonds, and o and • is with 9 x 10 6 bonds. 



6 = 2 self-dual 0.8915(6) 0.991(4) 

6 = 3 self-dual 0.8903(2) 0.998(1) 

6 = 4 self-dual 0.8892(6) 1.005(4) 

6 = 5 self-dual 0.8895(6) 1.003(4) 

6 = 6 self-dual 0.8890(6) 1.006(4) 

6 = 7 self-dual 0.8891(6) 1.005(4) 

6 = 8 self-dual 0.8889(6) 1.006(4) 



TABLE II: The locations of the multicritical points for the ± J 
Ising model on the self-dual hierarchical lattices. Also shown 
are the values 2H(p c ), which should be unity according to the 
conjecture. 



Lattice 



Pi 



P2 H(pi) + H{p 2 ) 



1 


0.9338(7) 0.8265(6) 


1.017(4) 


2 


0.8149(6) 0.9487(7) 


0.983(4) 


3 


0.7526(5) 0.9720(7) 


0.991(5) 


4 


0.8712(6) 0.9079(6) 


0.998(4) 


5 


0.8700(6) 0.9081(7) 


1.000(4) 


6 


0.9337(7) 0.8266(6) 


1.017(4) 


7 


0.9084(6) 0.8678(6) 


1.005(4) 


8 


0.9065(6) 0.8686(6) 


1.009(4) 



TABLE III: The locations of the multicritical points for the 
±J Ising model on mutually dual pairs of the hierarchical 
lattices. The results for lattices number 1 to 3 reproduce the 
results by Hinczewski and Berker. 



gations. In the thermodynamic limit, all the plots as in 
Fig. [6] become step functions. However, we investigated 
the finite-size hierarchical lattices. The slopes of all the 
plots are finite. We have checked these finite-size effects 
as shown in Fig. [7] From these analyses, the final error 
bars have been chosen to be p c /VNb for ± J Ising model 
and Jo/ 1 J 2 \/Nb for the Gaussian model, where Nb is the 
number of bonds of the hierarchical lattice. 



V. REPLICATED SYSTEMS 

There are slight differences between the results by 
the conjecture and the numerical data by the renor- 
malization group analysis for the multicritical points of 
quenched systems as shown in the previous section. We 
examine the conjecture for replicated systems on the 
self-dual hierarchical lattices in this section. Because 
the conjecture is based on the duality and the replica 



The results for the ± J model are given in Tables |TT] 
and IIII1 and those for the Gaussian model are in Table 
IIVI We show the values of the binary entropy H(p) for 
comparison with the conjecture. For the Gaussian Ising 
model, we also give the values of Hc{Jq/ J 2 ) for the self- 
dual hierarchical lattices. Similarly to the case of the ± J 
Ising model, it can be shown that the summation of both 
of values of H G { J / J 2 {A)) and H G {J /J 2 {B)) should be 
unity for the mutually dual pairs. We show such val- 
ues for comparison in Table IIVI We express here pairs 
of the locations of the multicritical points as Jo /J 2 (A) 
and J /J 2 (B). Seeing these results, we confirm slight 
but non-negligible deviations from unity for both cases 
of the ± J and the Gaussian models. We find the general 
tendency that the difference from unity for the Gaussian 
Ising model is smaller than in the ± J Ising model. 



Lattice Jo 2Hc(Jo) 



6 = 3 self-dual 1.0209(3) 1.0011(4) 



Lattice 


MA) 


MB) 


Hg(MA)) + Hg(MB)) 


1 


0.7605(5) 


1.3174(9) 


1.0005(8) 


2 


0.7655(5) 


1.3118(9) 


1.0000(8) 


3 


0.5569(4) 


1.6151(11) 


0.9999(7) 


4 


0.9704(7) 


1.0730(8) 


1.0009(10) 


5 


0.9701(7) 


1.0733(8) 


1.0009(10) 


6 


1.3175(9) 


0.7606(5) 


1.0003(8) 


7 


1.1450(8) 


0.9040(6) 


1.0001(9) 


8 


1.1436(8) 


0.9055(6) 


1.0005(9) 



TABLE IV: The locations of the multicritical points for the 
Gaussian model on the 6 = 3 self-dual and mutually dual 
pairs of the hierarchical lattices. 
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method, @, H, H, [HJ EH it is expected that we find such 
discrepancies also for replicated systems with the replica 
number n of natural numbers as in the quenched system 
(n-0). 

If the partition function is a single-variable function, 
we can obtain the transition point as the fixed point of the 
duality. Then equation xq(K) = Xq(K) gives the exact 
transition point, where Xq(K) and Xq(K) are the original 
and dual principal Boltzmann factors. 0, H, H, GlGjJ We 
illustrate this point by the pure Ising model as 

,K 



x (K) = 
xt{K) = 



V2 



(11) 
(12) 



By equating xq(K) and Xq(K), we find the transition 
point e~ 2Kc — \[2 — 1 for the pure Ising model on the 
self-dual square lattice. We assume that the equation 
xq (K) = Xq (K) is also satisfied at the multicritical point 
for the replicated systems as well as for the quenched sys- 
tem [n — ► 0), though the replicated systems have com- 
plicated interactions. 0, H H M [H 

Let us set K = K p , which is the condition of the Nishi- 
mori line, for the replicated ± J Ising model. The quan- 
tity K p is defined as e~ 2Kp = (1 — p)/p. Both of the 
principal Boltzmann factors are then given as,0, Q 

x (K) = 2cosh{(n + 1)K} (13) 
x*(K) = 2% cosh™ K. (14) 

Equation x$(K) = Xq(K) gives the conjecture ([2]) in the 
limit n — > 0. Validity of this conjecture can be rigorously 
shown for n = 1 and 2, and for n — 3, the same has 
been numerically confirmed for the square lattice. 3 It is 
worthwhile to examine whether Xq(K) = Xq(K) is sat- 
isfied or not at the multicritical point for the replicated 
± J Ising model on the self-dual hierarchical lattices. Be- 
cause the replicated ±J Ising model does not have any 
randomness for couplings, which has been taken into con- 
sideration by the configurational average, we can derive 
directly the multicritical points, by evaluating recursion 
relations such as in Eqs. ([5|) and ([6]). 

We obtained the locations of the multicritical points 
on the several self-dual hierarchical lattices from 6 = 2 
to b = 6 with the replica number n = 1, 2, 3, and 4. 
The results are shown and compared with the predic- 
tions by the conjecture x (K) — Xq(K) in Table IVl The 
results obtained in the previous section for the quenched 
system (n — > 0) are also shown for comparison. Equa- 
tion xq{K) — Xq(K) gives the exact answer for n = 1 
and 2. For n — 3 and 4, the multicritical point locates 
slightly away from the results of the conjecture. There- 
fore the assumption of the validity for the conjecture is 
violated for the self-dual hierarchical lattices. Consid- 
ering b — > oo, we find that the conjecture does not al- 
ways work well on this self-dual hierarchical lattice, since 
the system becomes an effectively one-dimensional chain 
without finite-temperature transition in the limit b — > oo 
although the conjecture gives the results independent of 
b. 



b 


n 


pc ^numerical 


Pc ^numerical 


2 


n — ► 


0.889972 0.8915(6) 


-0.0015(6) 




1 


0.821797 0.821797 







2 


0.788675 0.788675 







3 


0.769563 0.768851 


0.000713 




4 


0.757348 0.755451 


0.001897 


3 


n — > 


0.889972 0.8903(2) 


-0.0003(2) 




1 


0.821797 0.821797 







2 


0.788675 0.788675 







3 


0.769563 0.769022 


0.000542 




4 


0.757348 0.755942 


0.001406 


4 


n — > 


0.889972 0.8892(6) 


0.0007(6) 




1 


0.821797 0.821797 







2 


0.788675 0.788675 







3 


0.769563 0.769649 


-0.000086 




4 


0.757348 0.757763 


-0.000415 


5 


n — > 


0.889972 0.8895(6) 


0.0004(6) 




1 


0.821797 0.821797 







2 


0.788675 0.788675 







3 


0.769563 0.7705020 


-0.000939 




4 


0.757348 0.7601328 


-0.002785 


6 


n — > 


0.889972 0.8890(6) 


0.0010(6) 




1 


0.821797 0.821797 







2 


0.788675 0.788675 







3 


0.769563 0.771376 


-0.001813 




4 


0.757348 0.762313 


-0.004965 



TABLE V: Differences between p c by the conjecture equation 
Xq(K) — Xq(K), and p n umcricai by the exact renormalization 
analysis for the n-replicated ±J Ising model on several self- 
dual hierarchical lattices. For n — > 0, Pnumerical denotes the 
results obtained by Nobre's technique. 



VI. IMPROVEMENT OF THE CONJECTURE 

In this section we propose an improvement of the con- 
jecture, which reduces discrepancies observed in Table IVl 
We here consider the partition function for the replicated 
±J Ising model and its dual one and discuss their rela- 
tionship following Ref. [9]. The partition function for 
the replicated ± J Ising model on the Nishimori line is a 
multi- variable function of the Boltzmann factors as 

Z(K) = x (K) Nn z(u 1 ,u 2 , ■ ■ ■ ,u n ), (15) 

where the u r are the relative Boltzmann factors defined 
as 

x r {K) cosh{(n+ 1 - 2r)K] 



U r (K) = 



x (K) 



cosh{(n + 1)K} 



(16) 



Here r denotes the number of antiparallel pairs among 
the n pairs. The duality gives the following relationship 
between the original and dual partition functions, H 



x (K) Nb z(u 1 ,u 2 , ■■■ ,i 



(17) 
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FIG. 8: A schematic picture to consider the renormalization 
flow and the duality for the replicated ±J Ising model. 



where the u* 
fined as 



are the dual relative Boltzmann factors de- 



~{K) = 



{tanh r K (r — even) 
tanh r+1 AT (r = odd). 



(18) 



Figure [5] shows the relationship between the curves 
(ui(K), U2(K) ),•••, u n {K)) (the thin curve going 
through p c ) and (u^ (K) , u 2 (AT) ) , • • • , u* n (K) ) (the dashed 
line). For convenience, we show the projections on the 
two-dimensional plane (ui,U2). For T — ► (K — > oo), 
u r (K) — ► for any r, corresponding to the point F in 
Fig. [5] As T changes from to oo, the point represent- 
ing (ui(K), v,2(K)), ■ ■ ■ ,u n (K)) moves toward the point 
P along the thin line in Fig. [5J Then the corresponding 
dual point {u\{K), ul(K)), • • • , u* (AT)) moves along the 
dashed line in the opposite direction from P to F. 

If we were to consider a model with a single-variable 
partition function, the thin curve would overlap the 
dashed line, which fact would be reflected in the relation 
u*(K) — u r (K*). Solving this relation, we obtain the 
duality relation for the coupling constant K*(K). For 
example, the pure Ising model has the reduced Boltz- 
mann factors as U\{K) = cT 2K and u*(K) = tanhA". 
We obtain the duality relation eT 2K = tanhA" from 
u*(K) = U\(K*). However the replicated ±J Ising 
model is given by the multi-variable partition func- 
tion of (ui(K),U2(K)), ■ ■ ■ ,u n (K)). The thin curve 
(ui(K) , U2(K)) , ■ ■ ■ ,u n (K)) does not coincide with the 
dashed curve (u*(K), u^K)), ■ ■ ■ ,m*(AT)). Therefore we 
cannot identify the critical point with the fixed point of 
duality. 

A new point of view from renormalization group helps 
us proceed further. Let us notice two facts concerning 
renormalization group transformation, (i) The critical 
point is attracted toward the unstable fixed pint, (ii) The 
partition function does not change its functional form by 
renormalization for hierarchical lattice; only the values 



of arguments change. Therefore the renormalized sys- 
tem also has a representative point in the same space 
(ui(K), u 2 (AT)), • • • , u n (K)) as in Fig. [5J with the renor- 
malization flow following the arrows emanating from p c 
and d c to C, ph and di to P, and pi and dh to F. We ex- 
press such a development of relative Boltzmann factors 
at each renormalization step on the n-dimensional hy- 

perspace as (i4 , u 2 , ■ ■ ■ ,uS, ), where the superscript 
means the number of renormalization steps. The renor- 
malization flow from the critical point p c reaches the 
fixed point C, (u[°°\ u 2 , • • • , it^°°' ) ). On the other hand, 
there is the point d c related with p c by the duality. We 
expect that the renormalization flow from this dual point 
d c also reaches the same fixed point C because p c and d c 
represent the same critical point due to Eq. (|17p . 

Considering the above property of the renormalization 
flow as well as the duality, we find that the duality re- 
lates two trajectories of the renormalization flow from p c 
and from d c , tracing the renormalization flows at each 
renormalization. In other words, after a sufficient num- 
ber of renormalization steps, the thin curve representing 
the original system and the dashed curve for the dual 
system both approach the common renormalized system 
depicted as the bold line in Fig. [8j which goes through 
the fixed point C. 

It is expected that the partition function therefore be- 
comes a single-variable function described by the bold 
curve. This fact enables us to find the duality relation 
and identify the multicritical point by the following equa- 
tion, 



„(°°) 



(K) 



► (oo) 



(19) 



similarly to xq(K) — Xq(K) for the pure Ising model. 
We therefore have to evaluate Eq. (fP9|) , not the relation 
xq(K) — Xq(K) for the unrenormalized, bare quantities, 
to obtain the precise location of the multicritical point 
on the hierarchical lattices. Equation (fll?]) is expected to 
predict the exact location of the multicritical point for 
the hierarchical lattice. 

It should be noticed that the relation xq(K) = Xq(K) 
predicts values very close to numerical estimates in many 
cases of regular lattices as indicated in Table [TJ This 
means that the effects of renormalization are not large 
for those systems. If we regard the relation xq(K) — 
Xq(K) as the zeroth approximation for the location of 
the multicritical point, it is expected that the relation 



(K) 



(AT) is the first approximation and leads 



to more precise results than the relation xq(K) — Xq(K) 
does. We therefore propose the first-approximation equa- 



tion Xq 1 '' = Xq^ as the improved conjecture. We eval- 
uate the performance of this approximation in the next 
section for hierarchical lattices. 
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VII. RESULTS BY THE IMPROVED 
CONJECTURE 



In this section, we report the results by the improved 
conjecture Xq (K) = Xq^(K) and evaluate its perfor- 
mance compared with the conventional conjecture. 

The relation xo(K) = Xq{K) of the conventional con- 
jecture yields an equation that the binary entropy H(p) 
equals to 1/2 for self-dual hierarchical lattices as in Eq. 
(|2"|). Similarly to this relation, the improved conjecture 



(K) 



(K) gives an equation in terms of the 



binary entropy given by the values of the renormalizcd 
couplings as described below. After one-step renormal- 
ization, we obtain again the replicated Ising model on 
the hierarchical lattice with the renormalizcd couplings 
{K^} and their distribution function pW(ify). Here 
the renormalized quantities are determined by the initial 
condition. The original and dual principal Boltzmann 
factors for the replicated Ising model after one-step renor- 
malization are given as 



_|_ e -K i3 

71 



(20) 



(21) 



where the distribution function is given by, with the cou- 
plings {K^} obtained by such Eqs. (JSJ and ©, 

lit J 

(22) 

The product runs over the bonds on the unit cluster of 
the hierarchical lattice. Using these principal Boltzmann 
factors, we take the leading term of the replica number 
ri — * of the equation x^{K) = Xq^\K) and obtain 
the improved conjecture for the quenched system as 

J dKijPWiKij) log 2 {1 + cxp {-2Kij)} = ~. (23) 

The left-hand side of this equation will be written as 
HW(p). Equation (23]) gives the results for the replica 
number n — ► shown in Table IVII Similarly, we can ob- 
tain an equation to predict the multicritical point for the 
replicated ± J Ising model with a finite replica number n, 
whose solutions for n = 1, 2, 3, and 4 also shown in Table 
IVII All results are in excellent agreement with the nu- 
merical ones within their error bars. Comparison of Table 
IVII with the Tabic IVl clearly indicates significant improve- 
ments. We also find the improved conjecture gives results 
depending on the feature of each hierarchical lattice be- 
cause the prediction for the self-dual hierarchical lattice 



b Tl Pc pnumerical 


Pc ^numerical 


2 0.892025 0.8915(6) 


-0.0005(6) 


1 0.821797 0.821797 





2 0.788675 0.788675 





3 0.769048 0.768851 


0.000197 


4 0.755986 0.755451 


0.000535 


3 0.890340 0.8903(2) 


0.0000(2) 


1 0.821797 0.821797 





2 0.788675 0.788675 





3 0.769138 0.769022 


0.000116 


4 0.756250 0.755942 


0.000308 


4 0.889204 0.8892(6) 


0.0000(6) 


1 0.821797 0.821797 





2 0.788675 0.788675 





3 0.769629 0.769649 


-0.000020 


4 0.757619 0.757763 


-0.000144 


5 0.889522 0.8895(6) 


0.0000(6) 


1 0.821797 0.821797 





2 0.788675 0.788675 





3 0.769968 0.770502 


-0.000534 


4 0.758461 0.760133 


-0.001672 


6 0.889095 0.8890(6) 


0.0000(6) 


1 0.821797 0.821797 





2 0.788675 0.788675 





3 0.769947 0.771376 


-0.001429 


4 0.758300 0.762313 


-0.004013 



TABLE VI: The results by the improved conjecture Xq (K) 



is different from each other, which was not the case before 
as seen in Table IVl 

We can also see the performance of the improved con- 
jecture from another point of view. We can predict the 
phase boundary by the conventional conjecture if we do 
not restrict ourselves to the Nishimori line K p = K. The 
well-known transition point T c — 2.26919 for the pure 
Ising model is exactly reproduced by the conventional 
conjecture. However, except for this transition point, 
the conventional conjecture fails to derive the precise 
phase boundary especially below the Nishimori line as 
seen in Fig. [9l because the phase boundary of the ±J 
Ising model is expected to be vertical or slig htly reentrant 
below the multicritical point. 0, d, [ID, [24[ We find also 
inaccuracy of the conventional conjecture in the slope of 
the phase boundary at the transition point of the pure 
Ising model given as 



1 dT 

T c dp 



3.43294. 



(24) 



T=T C 



However this value is estimated as 3.23(3) by Nobre's 
method [l3| on the b — 3 self-dual hierarchical lattice 
and as 3.209 by the exact perturbation for the square 
lattice. [25| The improved conjecture, on the other hand, 
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0.88 0.9 0.92 0.94 0.96 0.98 



P 



FIG. 9: The phase boundary by the conventional and im- 
proved conjectures for a self-dual hierarchical lattice with 
6 = 3. The vertical axis is the temperature, and the hori- 
zontal axis is the probability for Jij = J > of the ± J Ising 
model. The bold dashed line is by the conventional conjec- 
ture and the bold solid line is by the improved conjecture. 
The thin dashed line is the Nishimori line. 



yields for b = 'i 



1 



dp 



= 3.30712. 



(25) 



T=T a 



This value is closer to 3.23(3). Also the phase boundary 
below the Nishimori line is modified as moving toward the 
p-axis as in Fig. [5J Thus the improved conjecture works 
better than the conventional conjecture for describing the 
phase boundary. 

The improved conjecture also succeeds in leading to the 
relation between the multicritical points on the mutually 
dual pairs. It is straightforward to apply the improved 
conjecture to the mutually dual pairs, similarly to the 
case of the conventional conjecture [l(| as, 



(26) 



where pi and p 2 denote the locations of the multicritical 
points on mutually dual pair. We estimate the values of 
the left-hand side of Eq. (|26|) for several pairs of hierar- 
chical lattices in Figs. [Handful The estimated results are 
given in Table IVIII We use the values of the locations 
of the multicritical points obtained by Nobre's method, 
as in Tabic IIII1 to compare the performance of the im- 
proved conjecture with that of the conventional conjec- 
ture H(pi) + H( P 2) = 1. There are cases in which the 



Lattice 


Pi 


P2 


value 


1 


0.9338(7) 


0.8265(6) 


1.002(7) 


2 


0.8149(6) 


0.9487(7) 


0.984(9) 


3 


0.7526(5) 


0.9720(7) 


0.993(9) 


4 


0.8712(6) 


0.9079(6) 


1.007(6) 


5 


0.8700(6) 


0.9081(7) 


1.011(6) 


6 


0.9337(7) 


0.8266(6) 


1.003(7) 


7 


0.9084(6) 


0.8678(6) 


0.996(6) 


8 


0.9065(6) 


0.8686(6) 


1.003(6) 



TABLE VII: The results by the improved conjecture for the 
mutually dual pairs. We estimate values of the left-hand side 
of Eq. (|26|l by the improved conjecture, shown on the right- 
most column of this Table. 



improved conjecture agrees with the numerical estimates 
for the mutually dual pairs, for the lattices of type 1, 3, 
6, 7, and 8 hierarchical lattices in Figs. Q] and [5] Un- 
fortunately, we find three cases for the lattices of type 2, 
4, and 5 in which the value of the left-hand side of the 
relation (|26|) is not unity within the error bars. However 
we find impressive improvements in Table. IVIII compared 
with the previous results in Table. IIIII 



VIII. DISCUSSIONS 

In the present paper, we first showed the existence 
of slight differences between the conventional conjecture 
and the numerical results for the locations of the mul- 
ticritical points on several hierarchical lattices. These 
discrepancies for the quenched system are caused by vio- 
lation of satisfaction of the equation xo {K) = Xq (K) for 
the replicated systems. This equation xq(K) = Xq(K) is 
satisfied for the case that the partition function is written 
by a single variable as in the pure Ising model. We expect 
that the partition function can be written as a single- 
variable function after a sufficient number of renormal- 
ization steps, considering the fact that the plot describ- 
ing the original model overlaps that of the dual model 
as in Fig. [8] Based on this consideration, we proposed 
the improved conjecture as the first approximation of the 
exact relation to determine the critical point. Through 
the derivation of the improved conjecture, one finds that 
the multicritical point on the self-dual lattice is given 
as a special point where the binary entropy given by the 
renormalized values on the unit of the hierarchical lattice 
becomes one half. If we need the very precise location of 
the multicritical point, we may use the numerical meth- 
ods for the renormalization group analysis to evaluate 
Eq. ((Til). 

The present study also gives a basis for the improve- 
ment of the conjecture on regular lattices. The improved 
conjecture for the hierarchical lattice reflects individual 
characteristics of each hierarchical lattice, because it in- 
cludes the renormalized couplings and corresponding dis- 
tribution function, which depend on the structure of the 
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hierarchical lattice under consideration. Similarly to the 
case of such hierarchical lattices, if we adequately carry 
out the renormalization for regular lattices, it should be 
possible to improve the conjecture also for regular lat- 
tices. Work in this direction is in progress. 
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